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Abstract 

We consider an irreversible autocatalytic conversion reaction A + B ^ 2A under subdiffusion 
described by continuous time random walks. The reactants' transformations take place indepen- 
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^ ' dently on their motion and are described by constant rates. The analog of this reaction in the case 

of normal diffusion is described by the Fisher-Kolmogorov-Petrovskii-Piskunov (FKPP) equation 

H ! leading to the existence of a nonzero minimal front propagation velocity which is really attained 

'^ \ by the front in its stable motion. We show that for subdiffusion this minimal propagation velocity 
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is zero, which suggests propagation failure. 
PACS numbers: 05.40.Fb, 82.40.-g 



I. INTRODUCTION 

The theory of reactions under subdiffusion had attracted recently considerable attention 
both because of theoretical and mathematical challenges posed by such problems, and also 
due to their growing practical relevance for description of phenomena taking place in porous 
media(as exemplified by geophysical structures) and in crowded interior of living cells (e.g. 
Refs. [1-3]). Several recent works were dedicated to the theoretical description of Turing 
patterns and of fronts in such systems [4-7]. Thus Ref. [7] concentrates on the front behavior 
in the system which in the case of the normal diffusion would be described by the Fisher- 
Kolmogorov-Petrovskii-Piskounov (FKPP) equation and shows, that under the description 
adopted, there exists a (minimal) stable propagation velocity of such front, just like it is the 
fact under normal diffusion. As we proceed to show, this is not always the case. Considering 
the fully irreversible analog of the reaction discussed in Ref. [7] under conserved overall 
concentration we show that the minimal propagation velocity is zero, which corresponds 
essentially to propagation failure. The system we consider here corresponds to the A + B — >■ 
2A irreversible reaction, which in the case of the normal diffusion is also described by the 
FKPP equation. 

The FKPP equation [8], proposed by R. A. Fisher in 1937 [9] as a model for propagation 
of a favorable gene in a population, corresponds to a mathematical description of the (irre- 
versible or reversible) reaction whose main stage is a bimolecular autocatalytic conversion 
A + B ^ 2A. Initially the whole system consists of particles (individuals) of type B. The 
introduction of the A-individuals into some bounded spacial domain (which is described by 
an initial condition sharply concentrated in vicinity of the origin of coordinates) leads to 
a propagation of a front of A into the B-domain. Physically this corresponds to a front 
propagating into the unstable state. Under normal diffusion, the process is described by a 
partial differential equation 

dt dx"^ 

for an irreversible reaction, where the initial concentration of B is assumed to be homoge- 
neous and equal to B{x, 0) = Bq everywhere except for the vicinity of the origin. In this case 
the overall concentration is conserved, due to the local stoichiometry of the reaction which 
does not change the number of particles. Here and in what follows we denote the nature of 
particles in reaction equations by Roman letters, while the corresponding concentrations or 



particle numbers are written in italic. 

Using the conservation law, the corresponding reaction-diffusion equation can be rewrit- 
ten as 

The equation for the reversible reaction has the same form but different coefficients in front 
of the two last terms in the right hand side (e.g. [10]). The FKPP equation is the simplest 
model of front propagation into an unstable state and serves as a paradigmatic model for 
many related phenomena. 

The front velocity in the genuine FKPP equation is determined by that of its leading 
edge, i.e. by the behavior of concentrations for x -^ oo. Since the concentration of converted 
particles in this leading edge is very small, the equation can be linearized, and the possible 
velocity of the front is given by the analysis of the linear propagation problem [11, 12]. For 
the FKPP equation, Eq.(l) leads for x ^ oo (and B -^ Bq) to 

(or by the corresponding equation for 6B = Bq — B = A) having the exponential propagating 
solution A ~ exp(— a; + vt) for all v > 2\/ DkB^ where the condition on v is imposed by 
the natural boundary condition A = for a; ^ oo and by forbidding oscillatory behavior. 



Further analysis shows that the minimal velocity v = 2\uJkBQ is exactly the one attained, 
which fact is known as the marginal stability principle. The front in the FKPP system is 
an example of the so-called "pulled front" , as it is "pulled" into the unstable state by its 
leading edge, and its propagation velocity does not depend on what happens in the interior 
of the front where the conversion of the most particles takes place. 

Modeling reaction phenomena in subdiffusive media requires some preliminary assump- 
tions on the nature of the transport process. Parallel to Refs.[13, 14] we assume that the 
subdiffusive motion on a mesoscopic scale is a consequence of trapping of particles on a meso- 
scopic scale, due to e.g. the bottlenecks connecting the voids in a porous system, while on a 
microscopic scale within the pores the A + B ^ 2A reaction takes place in a homogeneous 
solution and follows the mass action law. The local conservation of particle concentration is 
an inherent property of such systems which also holds on the mesoscopic scale. The model 
adopted corresponds therefore to the systems which at smaller scales consists of compart- 
ments in which the reaction follows the usual kinetic laws, while the subdiffusive transport 



between the compartments is described by continuous time random walks with a distribu- 
tion of sojourn times which may lack the first moment. As we proceed to show, due to 
the coupling between the reaction and transport term inherent for subdiffusion [13, 14], the 
behavior of the reaction front under subdiffusion and under the conditions discussed above 
is vastly different from the one under normal diffusion. The minimal propagation velocity of 
the front is zero, which corresponds to propagation failure. The preliminary results of our 
numerical simulations confirm that the front's velocity decays with time. These numerical 
results will be discussed in detail elsewhere. 

II. THE A + B ^ 2A REACTION UNDER SUBDIFFUSION 

A. General considerations 

Following the same procedure as in [4, 13] one obtains the equations for the particle 
concentrations as a consequence of the conservation laws. Discretizing the system into 
compartments (sites) numbered by the index i we get for the mean numbers of B-particles 
at site i 

Ht) = ^j7-i W + Ijr+iit) - jr(t) - kAit)B,it). (3) 

Eq. 3 is a local balance equation for the number of B-particles at site i with j~{t) being loss 
fiuxes of particles B from the site i at time t given by 
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The first term on the right hand side of Eq.(4) for the corresponding fiux gives the con- 
tribution of those B-particles that were at i from the very beginning and survived until t. 
The second term describes the particles, that arrived at i at a time t' and did neither react 
nor perform a jump until t. Equation (5) gives the survival probability of B-particles. The 



A-concentration is then given by 

Mt) 
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fr{t) = mPAit,0)B,iO)+mAiO) 

+ f ^p{t- f) U(t') + /r(t') - kA{t')B,{t') 



(6) 



dt' 



+ ^{t-t')PA{t,t') B,{t')+jr{t') + kA{t')B,{t') 



dt' 



PA{t,t') = I-PbHX) 



(7) 
(8) 



The loss fluxes for A-particles are denoted through fi{t). The first term on the right hand 
side of Eq.(7) corresponds to the particles that were at i from the beginning and converted 
from B to A until t. The second term represents the A-particles that were at i from the very 
beginning. The third and fourth term describe the particles that arrived at i at a time t' as 
A-particles, or arrived as B-particles and reacted until t. The probability Pa to gain new 
A-particles arises from the conservation of the total number of particles and the probability 
Pb that the B-particles react. The A-concentration depends on the B-concentrations at all 
previous times. 

Equations (3-8) are consistent with the conservation of the total number of particles in 
the reaction A + B ^ 2A. Let Ci = Ai + Bi be the total mean particle number at i and 
g^{t) = fj^{t) + jj^{t) the total loss flux of particles at i. From Eqs. (3, 6) we obtain the 
balance equation: 



a(t) = gtit) - giit) = l(?r i(t) + ]^gr^^{t) - gr{t). 



(9) 



From Eqs. (4, 7) follows that 

gr{t) = i,{t)[Bm + Ar 
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This equation can be solved by means of the Laplace transform, 

utp{u) 



dt'. 
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1 - ^(m) 

which yields a diffusion equation for Ci(t) in the Markovian case i/j(t) = ^ exp [—7], and a 
subdiffusion equation in the non-Markovian case iplt) oc t~^~", i.e. shows that the behavior 



of the total particle concentration is diffusive or subdiffusive, respectively. Moreover, if we 
choose the initial condition in a way that Cj(t) = const ^ the total number of particles is also 
locally conserved, Ci{t) = Ci{0) = Cq. This reduces the overall problem to the one for only 
one species. Passing to the continuous variables x = ai we obtain: 

B{x,t) = —Ar{x,t) - k{Bo - B{x,t))B{x,t) 
j-{x,t) = tp{t)PB{x,t,0)B{x,0) + / ^{t-t')PB{x,t,t') X 



B{x, t') + j-{x, t') + k{BQ - B{x, t'))B{x, t') 



dt' 



PB{x,t,t') = exp 



-k / {Bo-B{x,t"))dt" 
If 



(12) 



where the concentration of B is given by B{x,t) = Bi{t)/a^, and the corresponding dimen- 
sional constant is absorbed into a new reaction rate k = na^. The concentration Bq is set 
to unity in what follows. Eq.(12) can be rewritten in a form following from the equations of 
Ref. [15] for the irreversible case: 



B{x,t) = -k{l - B{x,t))B{x,t) + —A / M{t-t') 

2 Jo 



xB{x,t') exp 
with M{u) = wp{u)/[l - ^p{u)]. 



- I k{l-B{x,t"))dt" 
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dt' 



(13) 



B. Leading edge linearization 



To analyze the behavior in the leading edge we note that A{x,t) = 1 — B{x,t) becomes 
small for a; ^ oo and hence 

^ = — f A!^M{t-t')B{x,t')exp -k I {l-B{x,t"))dt" \ dt' 

~k{l-B{x,t)) (14) 



where we have interchanged the sequence of differentiation over x and temporal integration. 



For a; — > oo we have exp 
the form 
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1 so that the integrand can be put into 
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We now assume B{x, t) to he 1 — Aq exp [—X{x — vt)] for large x, i.e. consider a wavefront 
with an exponentially decaying leading edge moving at a constant velocity v. Inserting this 
into (14) and retaining only the terms of the first order in Aq, we get 
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M{t - t') exp [-X{x - vt')] dt' 



AQexp[-X{x-vt)] / M{t-t')dt', (15) 
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with the memory kernel M{u) = uiIj{u)/{1 — ipiu)). We note that the traditional way to 
proceed is first to linearize the equation and then to insert the exponential solution, however, 
since the solution of the linearized equation is anyhow an exponential we can also proceed 
the other way around, saving on tedious calculations. 

a. Markovian case Let us first show that for the Markovian case of exponential waiting 
time distribution the standard expression for the minimal velocity of the stable propagation 
is reproduced. Taking ip{t) = ^ exp [— ^] , one obtains M(t — t') = ^6{t — t') and 
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= \v exp [—\{x — vt)] — kAf) exp [— A(a; — vt)] \^Aq exp [— A(a; — vt)] 

It 

which ioY z = X — vt^ z —^ oo leads us to the following equation: 
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Y -v\^k = ^. 



This is the standard dispersion relation for the FKPP front, showing that real values of 
A (corresponding to physically sound solutions with nonnegative concentrations) are only 
possible for v > Vmi-. 



= 1^ d^kjlT = 2\/ Dk with D = a^/2r being the diffusion coefficient 
(note that Bq is set to unity). Note that the corresponding result emerges due to the 
cancellation of two terms of different nature in Eq.(16), which, as we proceed to show, does 
not take place in the non-Markovian case. 



b. Non-Markovian case We now assume the waiting time pdf to decay as a power law 
^{t) oc t"^"", < a < 1 for large t. With i = t - t' we get 

/ M{t - t') exp [\vt'] dt' = exp [Xvt] / M{i) exp [-Xvi] di 
Jo Jo 

= exp[\vt]M{\v). 



and moreover that 



M{t-t') exp [Xvt] dt' 



const 



-t 



a-l 



SO that the last term in (15) vanishes for large t. Altogether we have then: 
d 



dt 



(l — Aq exp [— A(a; — ft)] ) = —kAoexp[—\{x — vt)] 
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and with z = x — vt 
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For example, taking ilj{t) to be given by a Mittag-Leffler function ■?/'(t) = i5Q,[— (t/r)"], we 
have exactly ip{u) = [1 + (-ur)"]^^ and M{u) = t~°'u^^°', so that 







This equation is equivalent to 
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(18) 



and always possesses two nonnegative roots. Therefore the minimal propagation velocity in 
this case is zero. The propagation failure corresponds essentially to a continuum approaching 
of Vmin oc vD to zcro for the case when D ^ as it is the case in subdiffusion. 



III. COMPARISON WITH A RELATED MODEL OF REF.[7] 

Our result here differs from the one of Yadav et al.[7], where the stable front propagation 
with the constant velocity was found. The authors of Ref . [7] consider a situation in which 
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the A-particles undergo a reversible reaction corresponding to the branching-coalescence 
scheme, 

A + B -^ 2A 
2A ^ B + A, 

which for the Markovian case is also described by the FKPP equation. The equation used 
in Ref. [7] for the non-Markovian case reads: 

dA{x,t) 



kA{x,t) -ki[A{x,t)y 



dt 

+ a^A\ I M[t-t']A{x,t') xexp 



dt' \ . (19) 



kA{x,t")dt'' 
'0 L Jt' 

The velocity of the front was then obtained by means of the Hamilton- Jacobi approach via 

a hyperbolic scaling (a more elegant method known to lead to the same results as leading 

edge linearization in the Markovian case) and reads in our notation 




k{a — 3) 

r 



a 



"°^=r"°^^'"(^)"°<2-")- p"' 



with a generalized diffusion constant Ka = a^/r"'. 

Since the term quadratic in A does not contribute to the linearized solution determining 
the pulled front properties, the only difference between our approach and the one of the 
authors of Ref. [7] is that Eq.(19) is written for the concentration of the A-particles which 
are created during the reaction, while our equation Eq.(13) is put down for the B-particles 
which irreversibly disappear. Taking Eq.(19) for granted, we first checked that the method 
of the leading edge linearization used in the previous sections yields the same result for the 
minimal front velocity as in Ref. [7]. Therefore the differences between the results are not 
connected to the method applied, but to the equation itself. 

The difference in the results is explained by the fact that the equations for particles which 
are created during the reaction under the conditions discussed look differently from Eq.(19): 
Since the particle making its step as A might have entered the corresponding site both as 
A and as B (and having undergone a transformation), the equation for the product of the 
reaction (in our case A) has to contain the Laplacians of both A and 5-concentrations, as 
discussed in [13]. Moreover, considering reversible reactions put additional problems due to 
possible multiple A ^ B transformations during the waiting time at one site [14] . Although 
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leading to the classical FKPP-results in the Markovian case, the equations proposed in 
[15] are not well suited for the description of particles created during the reaction and for 
reversible reactions under conditions discussed in the beginning of the paper. 

IV. CONCLUSION 

Let us summarize our findings. We considered the reaction-subdiffusion problem for the 
irreversible autocatalytic conversion reaction A + B ^ 2A which, in the case of normal 
diffusion, is described by the FKPP equation. We show that in the case of subdiffusion, 
under the assumption of the reaction taking place irrespectively on the particles' mesoscopic 
motion, the minimal front's propagation velocity is zero, suggesting the propagation failure. 
This regime can be considered as a result of a continuous transition to -D ^ in the normal 
FKPP situation. The result obtained is in contrast to the treatment in [7] where the finite 
propagation velocity of a traveling wave was found. The differences in the assumptions of 
the two approaches (corresponding essentially to two very different situations) are discussed. 
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